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A systematic search for possible order parameters of the so-called hidden order of URu2Si2 is 
conducted. Among the possible candidates that do fulfill the experimental symmetry restrictions 
on the hidden order parameter, we find one candidate that stand out - one of the components of 
a triakontadipole multipole tensor that belongs to the Ai u irreducible representation of the point 
group Dih- This solution is characterized by a Q = (0,0, 1) ordering of the triakontadipoles and 
has a symmetry that forbids magnetic moments as well as most other multipole ordering on the 
uranium sites. This hidden order phase is closely related to the the antiferromagnetic phase, which 
is manifested in the similarities of the geometries of the calculated Fermi surfaces. Finally is is 
found that this non-magnetic solution allows for another secondary superimposed order parameter 
that belong to B^u, which gives rise to an anisotropic in-plane susceptibility without breaking the 
tetragonal crystal symmetry. 

The hidden order (HO) of the heavy fermion material URii2Si2 has attracted a lot of interest since its discovery 
in 1985 p], as described in a recent review [2]. The main enigma is at 17 K there is clear signature of a secondary 
phase transition but not of any observable order parameter (OP) in the HO phase. The only signal is that there is 
a tiny staggered magnetic moment which cannot account for the change in entropy at the phase transition. Another 
important aspect of the HO is that it under pressure goes through a phase transition to an anti-ferromagnetic (AF) 
phase with relatively large moments. This phase transition is now conclusively established to be of first order [3J. 
This and the fact that the order parameter of the HO phase does not cause any change in observable symmetries, 
such as crystal symmetry or magnetic moments, put some severe conditions on its symmetry properties. 

There are numerous theoretical suggestions for the HO parameter in the literature, see Ref. [2j for a fuller account 
of these. Some are explicitly itinerant in nature as e.g. unconventional density waves [I], orbital currents [5], or 
helicity order [5J. Other focus more on the local atomic order parameter. In fact all multipolar order up to rank five 
have been suggested; quadrupoles @], octupoles [5], hexadecapoles and triakontadipoles [TU] . 

In this Letter we will allow for a general real space atomic order parameter that can arise due to ordering of the 
itinerant uranium /-states. First we systematically study what symmetries of the OP that are consistent with various 
experimental observations. Then we test which of the OP candidates can be stabilized in a realistic calculations of 
the correlated /-electrons by means of a combined density functional theory and correlation treatment (DFT+J7). 
These complementary studies point towards which OP are allowed by symmetry and compatible with the electronic 
structure of URu2Si2- This study point conclusively towards an OP which is a superposition of two triakontadipole 
components belonging to different irreducible representations of the isogonal point group: A\ u B2 U - 

We will allow for a staggered OP as a superposition of several independent OP, where each OP takes a general form 

1 N - - 

n 

where Q is an ordering wave vector, R n are the uranium atomic positions, N the number of atoms in the crystal, f\ 
is the /-electron creation operator at atom site n and I Q is the operator for the local multipole of type a. Here Y a 
is a matrix-operator in the 14-dimensional space of /-orbitals and f n is a vector-operator in the same space. 

In order to have Eq. Q to include all possible OP stemming from the /-shell, a should enumerate all possible 
degree of freedom within this shell. This is known to be handled by the so-called tesseral multipole tensor moments 

(tmtm) umaiBnnina, with a = {k Pr -t}, 

w k r{n) = Trr^ </„/t) = (/t r**7 n ) , (2) 

where the corresponding expansion matrices T^ pr are known matrices in the /-orbital space, see also Supplementary 
Materials (SM) [TBJ for more details. The multipole tensors in Eq. ^ have a simple physical interpretation; for even 
k, they are multipoles of the charge (p—0) or spin-magnetization (p=l), while for odd k they are multipoles of the 
corresponding currents. The rank of the tensor is given by r and its time reversal (TR) symmetry is given by (— ) fe +P. 
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TABLE I: The TMTM components that corresponds to IR of the group Dah- In the enumeration of the components t, n > is 
an integer. The experimental compatibility for HO OP belonging to different IR of D^h are ranked from the discussion points 
i-iv. Plus (minus) sign means that it has some (dis-) advantageous features for TR-even (g) and TR-odd (u) IR, respectively. 
For point in there is the possibility of superimposed OP and the different letters indicate without any ranking order which two 
IR have to contribute together. 
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It is easy to show that TMTM have simple transformations rules under point group operations. For a rotation 
with an angle (f> around its quantization axis c and for a rotation by 7r around the first perpendicular direction a they 
behave as 

lZ(c, 4>) Wf Pr (n) = cos(t(j)) Wt pr (n) + sm(t(f>) wj%"(n) 

TZ(a, tt) w^ r {n) = {-) t+r sgn(t) , (3) 

respectively. In Table [I] the irreducible representations (IR) of the TMTM are determined for the tetragonal crystal 
point group through the characters of its two generators, C4 = 7\L(c, it/2) and c 2 = lZ(a, ir), where a and c denote 
the lattice directions. 

Before we present our calculations we will reinvestigate which of the TMTM OP that are compatible with the most 
important constraints put on the HO OP that have been gathered by the huge amount of experimental studies. It is 
not possible to cover all experimental aspects 0, so we concentrate on experimental observations that have simple 
and direct implication on the symmetry aspects of the HO. 

i) Non-broken lattice symmetry. At high temperatures URti2Si2 has the crystal symmetry of the space group 
139 (I4/mmm) with its isogonal point group D^h (4/mmm). It is established that there are no lattice distortions at 
the HO transition at 17 K, i.e. the HO phase also belong to the tetragonal crystal class. Since the uranium atoms are 
situated at maximally symmetric Wyckoff sites the local site symmetry is also D^h- Then from Table [i] it is clear that 
if the local four-fold rotation remains, the OP has to belong to the A\ or A2 IR of the point group. However when 
the crystal symmetries are taken into account one can see that the broken local four-fold symmetry can be taken care 
of by crystal symmetry operations, which has been noted earlier for the case of quadrupoles 18 . The remedy is a 
four- fold screw axis generated by the non-symmorphic symmetry operation (c4||i|). In the cases of B\ and B 2 IR 
the tetragonal symmetry is recovered for an ordering wave vector Q = (001), while for E the corresponding wave 
vector is (00|). 

ii) Vanishing magnetic moments. Since it is now established that the phase transition under pressure from 
the HO phase to the AF phase is of first order [3] , the HO OP cannot belong to A 2u since it is the IR of magnetic 
moments along the c-axis. If we are looking for a symmetry reason for the non-existence of magnetic moments in the 
HO phase, the IR of the HO OP cannot neither belong to E u , the IR of the in-plane magnetic moments. 

iii) Low symmetry in-plane susceptibility. Recently there have been reported that the HO phase has a broken 
four- fold symmetry [5 . An analysis has shown that this arises from the OP squared that interacts with the applied 
magnetic field squared, which leads to a non- vanishing B 2g IR of the magnetic susceptibility [51 UB] - In Ref- [B] it was 
shown that only the OP belonging to E has a B 2g IR in the direct product representation of its square, B 2g £ E v ®E Vl 
v = {g,u}. 

Another option is that OP is a linear combination of two independent OP of different IR. Then since [16] A\ V ®B2 V = 
A 2v <8> B lv = B 2g , we see that the co-existence of these type of OP would also lead to the observed variation of the 
magnetic susceptibility. Note that the two OP have to have the same TR symmetry v, i.e. even or odd. 

iv) Good hideout. From Table |I] we can observe that even when the primary OP would have a high rank which 
makes it hard to observe directly, it will leave traces in terms of induced multipoles of lower ranks. In that sense there 
are two type of OP that have the chance to be better hidden than others. They are OP with t = —4 and either even 
r > 4 belonging to IR A 2 or with odd r > 5 belonging to IR Ai. 

Summary of candidates. Our survey is summarized in Table [T] It is a linear combination of OP that best fulfill 
the criteria from point i-iv, either A\ u © B 2u or A 2g © Bi g , while the best pure OP belongs to either Ai u or A 2g . 
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Electronic structure calculations. In order to determine which of the TMTM are compatible with the electronic 
structure, we will perform a systematic survey in terms of realistic calculations. Care has to be taken to allow for 
TMTM OP solutions as of Eq. Q that belong to the different IR of the isogonal group D^h as listed in Table [I] In our 
approach we break the symmetry by inducing staggered u>f pr on the uranium sites and determine the largest possible 
symmetry group that is compatible with their existence. Then by iteration we determine if this starting assumption 
converges to a non-trivial solution. In principle it is possible that two cases of induced OP of the same IR do lead to 
two different solutions, but they may of course also converge to the same solution. 

The electronic structure is determined with the DFT+C/ approximation within the APW+Zo method as implemented 
in the ELK-code [HI 12] • The calculations were performed in a similar way as earlier described [3] [TO] , and is presented 
in more details in SM [TO]. 

For each possible TMTM component of each IR as given in Table [I] we have started a calculation with a large value 
of the corresponding multipole. All calculations enforce a staggering wave vector Q = (001) for tensors with rank 
r > 1. It was found that one solution exists for each TR-labcllcd IR and these results are collected in Table [TT] for the 
case of U = 1 eV. In order to quantify the importance of different tensor components we have utilized the concept 
of polarization 7r^ pr [7J, which is a normalization independent quantity that directly measures the importance of the 
different contributions to the polarization of the density matrix. It is proportional to the square of the components 
of the TMTM w^ pr and all components except kpr = 000 add up to a total polarization 7r tot , see SM [TO] for more 
details. 

In the case of TR-cvcn all IR converged to the trivial un-polarized case of Ai g with only the rotational invariant 
tensors, w^ 00 and iUq 10 , being non-zero. They correspond to expectation values of the /-occupation and the operator 
£ ■ s, respectively. The calculated large value of Wq 10 leads to an enhancement of the intrinsic spin-orbit coupling, 
which brings the solution in to a relativistic regime with predominantly j = 5/2 occupation. The resulting OP is not 
actually staggered. 

When the TR symmetry is broken we find in all cases that components of the triakontadipole wf 15 dominate. In 
addition we find smaller contributions from all TMTM components that are allowed by symmetry as given in Table 
[TJ that is both TR-odd as well as TR-even. However, only two TMTM have significant polarizations, Wq W and w^ 15 
for all cases except A 2u and E u . In all these cases the w^ 10 have similar polarizations as in the TR-even case. 

The case of A 2u is the most complicated. Here there is a strong competition for the exchange energy between 
the magnetic dipole moments along c-axis and the two allowed triakontadipole components. The magnetic dipoles 
have significant polarization for all three MP variants, spin Wq , orbital wj 01 , and spin dipolar u>q U , while the 
triakontadipole Wq 15 has a contribution of similar strength but w® 15 has the largest polarization. Although all 
initializations lead to the same solution, the convergence is often extremely slow which indicates that there is a flat 
energy-landscape. This was observed in the earlier study where we first observed the triakontadipoles in the A 2u phase 
[10] . For the two dimensional IR E u we have performed two calculations where the OP is oriented along the in-plane 
symmetry directions (100) and (110), respectively. Both the two solutions are stable although different, indicating a 
strong anisotropy. As anticipated both cases posses non-vanishing magnetic moments. The local magnetic moments 
are 1.1 and 0.8 /xr, respectively. No polarized solution could be found in the case of tetragonal E u , i.e. with the 
ordering vector Q/2. 

For the TR-odd case we have also looked for the possibility of superimposed OP. Only one new solution was found 
when allowing for various combinations of superimposed OP, but this is one of the HO OP candidates that were 
singled out fulfilling all the experimental constrains. It is Ai u (B B 2u , with the ip_f(Q) two orders of magnitude larger 
than ip2 15 {Q)- The product of these two OP interacts with the global magnetic field in the a6-plane and gives rise to 
a nonzero B 2g IR of the magnetic susceptibility [16] , which explains the in-plane response of the torque experiments 
0. 

In Fig. [I] the variation of the Fermi surface (FS) sheets are displayed, from the uncorrelated cases to the finite 
U = 0.5 cV and 1.0 eV cases for both the two solutions with largest polarizations, the A\ u and A 2u . At U = they 
reproduce the FS of Ref. [55]. A 2u corresponds to their AF solution, and the, in this case unpolarized, A lu corresponds 
to their paramagnetic. In this limit, we observe that the FS of the two solutions are radically different, while at finite 
U they become surprisingly similar. Hence if we would identify the non-magnetic solution with a A\ u OP as the 
HO phase and the magnetic solution A 2u as the AF phase, we see that the calculated FS are in excellent qualitative 
accordance with recent Shubnikov-de Haas measurements [27j that observe only minor changes in the phase transition 
from HO to AF phase under pressure. In the calculated case we can see very large resemblances in topology as well 
as sizes between the FS of the two solutions at both U = 0.5 eV and U — 1 eV. Since the FS for the two different 
[/-values do not even have the same topology, the FS geometry is very U dependent. So it is left for a future study 
to calculate more in detail the U variation of the FS and to compare more quantitatively with experiments. 
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TABLE II: The different sofutions for calculations with (7 = 1 eV. For each IR label, the total polarization, the dominating 
( n t Pr > 1) TMTM components, their values and their polarizations are given. For E u the case with largest polarization, which 
is with the two independent components along the plane diagonal (110), is presented. 
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FIG. 1: The Fermi surface of URu2Si2 as a function of parameter U for the case of OP belonging to the IR Ai u and A2 U , 
respectively. The green/purple, blue/yellow and red/light-blue colors on the different sides of the Fermi sheets refer to the 
electron/hole side. 
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In this study we have performed a detailed and systematic survey of possible HO OP in terms of ordering within 
the uranium /-bands. A picture arises from some common features of Tables [T| and \H\ First we observe there are no 
indication at all in the calculations for a TR-even HO OP. Secondly, it is clear that the AF solution at high pressure 
is the IR A 2u where the OP in the calculation is a superposition of magnetic vector OP and TR-odd rank 5 TMTM 
components, with the last ones dominating. Thirdly, this phase competes with an almost pure TMTM OP of the 
kind ip^}i(Q) which belong to the IR A\ u . These two solutions, A 2u and A lu , are among the ones that have the 
largest polarizations, as shown in Tabic [Tl| The HO phase A\ u only allows a few multipolar components with tesseral 
component t = —4, where all except become small. This leads to an optimally hidden OP. It is further found 
that there exist a solution where this A\ u OP is superimposed with a smaller B 2u triakontadipole OP, which would 
explain the recently observed anisotropic in-plane susceptibility [5]. 

From Table |H] we can directly observe that the various solutions for the different IR all have a large contribution 
from the rank 5 TMTM OP ip 615 . These results are in good accordance with our earlier observations that these 
trikontadipoles play a large role for URu2Si2 in particular |10j and for the time reversal symmetry breaking for 
moderate or strong spin-orbit coupling in general [7]. The latter is summarized as Katts' rules. After completion of 
this study we became aware of a recent calculational study [25J which also identify the importance of large rank 5 
multipolcs in the electronic structure of URu2Si2- In those random phase approximation (RPA) calculations only the 
j = 5/2 states are included (see 16] for a discussion) and it was concluded that the magnetic r = 5 E u state is the 
best candidate for the HO phase. 

The support from the Swedish Research Council (VR) is thankfully acknowledged. The calculations have been 
performed at the Swedish high performance centers HPC2N and NSC under grants provided by the Swedish National 
Infrastructure for Computing (SNIC). 
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ORDER PARAMETERS 



In this study our main focus is on order parameters (OP) of the general form 

1 N 

n 



N 

V ^ 



' ^e^^(n), (4) 



where Q is an ordering wave vector, R n are the uranium atomic positions, N the number of atoms in the crystal 
and /J is the /-electron creation operator at atom site n. Y kpr is the operator for the local tesseral multipole tensor 
moment (TMTM) component, 

w^(n) = TiT^Pn = Trlf r (fjt) = (fi f n ) , (5) 

where the trace is over the / orbitals. For a /-shell < fc < 6, < p < 1 and \k — p\ < r < k + p, which constitute 
26 different multipole tensors, of which 13 are time reversal (TR) even and 13 TR-odd, and the total number of 
tensor components are 196 = 14 x 14. This then accounts for the full freedom of the 14-dimensional density matrix 
p n = (f n fn)- In the 14-dimensional space of /-orbitals, Y kpr is a matrix-operator and f n is a vector-operator. In a 
{jnij }-representation of the /-states (jj-basis) we have 



r t fcp r 2 =^^(-) fe+p+r <| a s p \ii^; m2 (6) 




7^ 2 ; m2 =(-) h - mi r[ 31 T n ) . (7) 

\ —mi t rri2 

Here £ — 3, s — 1/2, the (■■•)- and {. . . }-symbols are the Wigner-3j and -9j, respectively, Nk pr i is a normalization 
factor and [a...b] — (2a+l)...(26+l). [HE] The operator Tbrings a spherical tensor, which was used in earlier studies 
[H |3] , to a tesseral form 

![a t + (-)*a_ t ] /V2 = V2Ua t t>0 
a t t = Q . (8) 

i[a t -(-ya_ t }/V2 = V2Za ltl t<0 

The tesseral form is convenient when considering rotational symmetries as in the present study. The TMTM in Eq. Q 
have a simple physical interpretation; for even fc, they are multipoles of the charge (p=0) or spin-magnetization (p=l), 
while for odd k they are multipoles of the corresponding currents. The rank of the tensor is given by r and its time 
reversal (TR) symmetry is given by (~) k+p . 

Hence all possible OP stemming from the /-shell is covered by a superposition of OP in terms of TMTM of Eq. [4] 
For instance the OP of an ordinary spin density wave is given by - 0° 11 (Q). This can be easily seen since T? 11 = at, the 
Pauli spin matrices in tesseral form, i.e. o\ = a x , cx_i = cr y and &q = a z , respectively. Another example is ip2 12 (Q)i 
which is one of the staggered quadrupoles models suggested in Ref. |U 



Symmetry properties 

It is easy to show that TMTM have simple transformations rules under point group operations. For a rotation with 
an angle (j) around its quantization axis (which we denote z) and for a rotation by 7r around the first perpendicular 
direction (which we denote x) they behave as 

lZ(z, (/)) w\ pr {n) = cos(t(j)) wl pr {n) + sin(i</>) wj^(n) 

K(x, tt) w k t pr (n) = (-) t+r sgn(i) w k t pr {n) , (9) 
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respectively. Hence the rotational properties are determined by the rank r and the component t only. In Table I of 
the main Letter (ML) the irreducible representations (IR) of the TMTM are determined for the isogonal point group 
D^h through the characters of its two generators, C4 = 1Z(c,tt/2) and c 2 = lZ(a,n). In addition the TMTM behave 
under a TR-operation as 

Qw k r{n) = (-) k+r w k r{n). (10) 

As the local site group for the uranium atoms is equal to the isogonal point group, Table ML-I describes the local 
symmetry for the different IR. Then it is clear that only A\ and A 2 are compatible with the fourfold rotational 
symmetry of the tetragonal space group. However, this symmetry is recovered also for the other IR if super-cells are 
allowed for. This comes from the non-symmorphic group element (c4|i|j), a fourfold screw operation that connects 
the corner (n even) and body-centered (n odd) uranium sites in the original bet structure. It maintains the tetragonal 
symmetry also for the other IR, as for Bi and B 2 we get 

^(0) = -*(1) (11) 



2 2 2. 
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^ 2 

111 



<±\\\\\ <2 r (0)=< 2 r (2) (12) 



while for E we get (q is an integer < q < (r — l)/2) 



) <r(o)=^(i) (is) 



'2 2 2; M w - -M* 
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<r(o) = -<r(2) (14) 



1 1 1 



C4I222) >f(0) = -^](3) (15) 



^l^^) 4 <f(0)=<f(4i. (K.i 



where for odd t = [q] = 2q + 1, w kpr and w;5 r span the two-dimensional IR. The operations leading to a change of 



sign of the TMTM, i.e. Eqs. (|ll[) and (14), correspond to ordering vectors Q — (001) and Q/2, respectively. 



Susceptibility 

As discussed in the ML there are no experimental signature that the tetragonal crystal symmetry is broken in the 
hidden order (HO) phase. However, recently in measurements of the in-plane susceptibility on single URu2Si2 crystals 
an anisotropic component was detected. [S] This was subsequently analyzed by Thalmeier and Thakimoto [6 . Here 
we will extend their analysis and show that the experiments can be explained by a non-vanishing Bi g contribution 
to the in-plane susceptibility. The experiments are based on measurements of the torque on the sample in a uniform 
external magnetic field H that are constrained to the ab plane of the tetragonal crystal. This torque is given by 

T = [i B V X HxH, (17) 

where V is the sample volume, \ t ne susceptibility tensor and H = (H a , Hi,, 0) = H (cos <j), sin <f>, 0). The in-plane part 
of the symmetric x tensor is decomposed into three IR 

TZep [ X \\] = A lg 8 Bi g 8 B 2g = TZep [ Xaa + Xbb] ® Kep [ Xaa - Xbb] ® Kep [ Xab ] . (18) 

Then the torque is given by 

f = fi B VH 2 c I ^(xaa - Xbb) sin 2<j) - X ab cos 20 j , (19) 
i.e. with two-fold symmetric contributions from the B\ g and B 2g IR of the susceptibility, respectively. 



TABLE III: The IR decomposition of the direct product of two IR: the factors given by the row and column and the direct 
product in their intersection, v labels the TR symmetry, either even g or odd u. 
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Furthermore, these susceptibilities will have contributions from the HO, *&{Q), through its interaction with the 
magnetic field in the Landau free energy expansion 

TH = 9l[^(Q) 2 }B lg [H 2 }B lg + 92^(Q) 2 ]B 2a [H 2 ]B 2g + ■■■ 

= 9i[^{Q) 2 ] Blg \{Hl - H 2 )+g 2 [*(Q) 2 ] B2g H a H b + ..., (20) 

where only second order interactions appear due to the staggering of the HO, Q 7^ 0, while the magnetic field is 
uniform. Thus the B\ g susceptibility are proportional to [5 , (Q) 2 ]s lg , i.e. the B\ g part of the IR decomposition of 
the squared OP, while the B 2g susceptibility are proportional to [^{Q) 2 ]B 2g , since by definition Xij — d 2 F /dHidHj. 



From Table III one can see that there are three ways each to get B\ g and B 2g contributions in the IR decomposition 
of a squared OP, i.e. 1Zep[*S>(Q)\ <X> TZep[^f(Q)]. B\ g can be obtained from A\ v ® B\ Vl A 2v ® B 2v or E v ® E v , while B 2g 
can be obtained from A\ v ® B 2u , A 2v ® B\ v or E v ® E Vl where v indicates the TR symmetry, even (g) or odd (u). 

In the torque measurements [5 a cos 2(f> oscillation was observed in the HO phase only, which imply that there is a 
non-vanishing \ab, i-e- a B 2g IR, in the presence of the HO OP. Hence, if this OP belongs to a single IR of the point 
group D ih , only TZep\^{Q)] — ^e,AQ) gives rise to the observed two-fold in-plane susceptibility Xab, i-e. only OP of 
IR E g or E u are candidates for the HO. On the other hand if the OP is a superposition of components belonging to 

different IR, OP of type TZep[^f{Q)] = * a 2 AQ) ® ^ b 1v {Q) ® ■ ■ ■ or ^ep[*(Q)] = 4%„(<3) © *s 2 „(Q) © ■ ■ ■ are also 
possible. Hence in a general case, with an OP having contributions from all possible IR, we have that 

d 2 T 

Xab = Q HadHb = J2 2 92 (*A lK *B aB + *A a „*B 11( + * E ua ^ E„ b ) ■ (21) 

Here *&E va and *&E ub span the two-dimensional IR E v and are chosen such that c^E va = E^ , c 2 E va = (—) r+l E ua 
and c 2 E ub = (—) r E va , with r the rank of the OP tensor. 



Relativistic effects 



In presence of strong spin-orbit coupling, there will be a large splitting between the j = £ — 1/2 = 5/2 and 
j = £+ 1/2 = 7/2 states. Then it is useful to study the density matrix in a jj-basis 

( 5 5 5 7 \ 

7 5 7 7 ) 

where each sub- matrix p^ n is spanned by m\ and m 2 with — ji < mj < j\ and —j 2 < m 2 < j 2 . 



Strong spin-orbit coupling 

Let us first discuss the limit of very strong spin-orbit coupling. In this limit the TMTM Wq 10 is related to w® 00 
through wj 10 = — |uio°°, and only the j = 5/2 occupation is non-zero. As we well discuss below, URu2Si2 does not 
fulfill this criterion and hence rather possess an intermediately strong spin-orbit coupling. However it is interesting 
to study this limit not least since it is assumed in some other theoretical studies. In this relativistic limit the j — 7/2 
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states are much higher in energy than the j = 5/2 states and the occupation is restricted to the sub-matrix p 
Eq. J22l. In this limit the TMTM are then given by 



kpr 

w t 



(n) ss TrT t pr p% ; 



N, 
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n Ol 



(22) 
(23) 



Now from the definition of 7j jr in Eq. (jjj, one can see that r in Eq. (23 1 is given by the vector coupling of two j = 5/2 
angular momenta, and hence can take values in between and 5. From general relations for exchange of two columns 
in the 9j symbols 




2l+2s+2j+k+p+r 




(24) 



Now the fact that 2(1 + s + j) is always even the 9j-symbol in Eq. ( 23 1 has to be zero for odd k+p + r. Hence only 
even k+p + r TMTM contribute in this limit. This in turn gives that odd (even) r multipolar moments have to be 
TR-odd (TR-even). 



The different rank r TMTM are then related through Eq. <[23j) , e.g. in the case of r = 5 the TMTM with kpr = 505, 
415 and 615 have fixed ratios of 
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0.013 



(25) 



-0.077. 



(26) 



In this case the 615 TMTM is largest which is just an example of the general case. The spin dependent (p = 1) 
TMTM have largest weight for the r = k — 1 tensor moments, as these are favored by the spin-orbit coupling for less 
than half-filled /-shell. 



Intermediately strong spin-orbit coupling 

For a full /-shell the TMTM OP with highest rank is the Vt 17 (Q) 01 rank 7 and in tm s case the TR symmetry is 
not given by rank r, as also odd k + p + r TMTM are allowed. They are included in all the analysis for completeness, 
although these tensor moments are somewhat obscure. Within these the orbital and spin degrees of freedom couple 
in an axial way, as of e.g. w 111 = %t x s. These tensor moments only arise from the off diagonal j\ 7^ j'2 blocks of 



Eq. (22 1 as they are not allowed in the block-diagonal part as the corresponding 9j-symbols have to vanish according 



to Eq. ( 24 1 . They will be included in the analysis and it would be fascinating if they would play a role, but in the 
present case of URri2Si2 they are of marginal interest. 

That URu 2 Si2 belongs to the case of intermediately strong spin-orbit coupling can be directly seen from the 
occupation numbers. Both the /-occupation, Wq 00 , as well as the effective spin orbit coupling, u>q 10 , are essentially 
independent on the assumed IR and take values around 2.6 and —2.5, respectively. This corresponds to that Wq 10 has 
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a value 70% of its saturation value, 



4, ,,ooo 

'3 W ■ 



The occupation number for the 5/2 sub-shell is then given by 



while for the 7/2 it is 
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000 
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(27) 



(28) 



The extra flexibility of having a non-zero 77.7/2 can lead to a stronger 615-polarization on the expense of the 415 
and 505 ones. This is e.g. confirmed in that the calculated ratios corresponding to Eqs. (25) and (26) are 0.006 and 
—0.052, respectively, for the case of A lu and t — —4. 



Polarization 

In order to quantify the importance of different tensor components we have utilized the concept of polarization [7] , 
■Kf W is a normalization independent quantity that directly measures the importance of the different contributions to 
the polarization of the density matrix. It is proportional to the square of the components of the TMTM w t pr and all 
components except kpr — 000 For each tensor component except, k = p = r = 0, 

= [eskpr]\N kpre w^ r \ 2 , (29) 

which all add up to a total polarization 7r tot = X^fcprt 71 *^ '' ^ ne total polarization is constrained by the inequality 
7r tot < nrih, where n is the occupation number of the /-shell and rih the corresponding number of holes. 



ELECTRONIC STRUCTURE CALCULATIONS 



In ML a systematic study which of the TMTM are compatible with the electronic structure were conducted through 
realistic calculations. 

The electronic structure was determined with the DFT+C/ approximation within the APW+/0 method as imple- 
mented in the ELK-code. [9] In the DFT+/7 approach a screened Hartree-Fock interaction is included among the 
5/ states only. The Slater parameters are determined individually by using a screened Yukawa potential, this 
approach has been showed to be particularly convenient since all four are controlled by the choice of a single 
parameter, the parameter U. [5] For double-counting we adopt the automatic interpolation scheme of Pethukov et 
al. [10] The calculations were performed at the experimental lattice constants. [T] The muffin-tin radius i?^T °f U 
is set to 1.9 A, and those of Ru and Si are set to 1.2 A. The parameter i?^ T |G -I- fc| max , governing the number of 
plane waves in the APW-Mo method, is chosen to be 9.5. To allow for the Q order the unit cell is doubled and the 
corresponding Brillouin zone is sampled with 18 x 18 x 10 fc-points. 

Care has to be taken to allow for TMTM OP solutions as of Eq. Q that belong to the different IR of the isogonal 
group Z?4/j as listed in Table I of ML. In our approach we break the symmetry by inducing staggered w^ pr on the 
uranium sites and determine the largest possible symmetry group that is compatible with their existence. Then by 
iteration we determine if this starting assumption converges to a non-trivial solution. 
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